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lution of  exponentially  distributed  random  variables.  After 
reducing  the  system  scenario  to  its  shorthand  format,  one 
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I-    INlSQDgCTIQW 


The  reliability    shorthand  considered    in    this    papsr   has 
been    developed    in   conjunction   with    the    course    Ok    4302    "Reli- 
ability   And   Weapon    Systems    Effectiveness    Measurement"    at   the 
Naval    Post-graduate    School.       A  tutorial    introduction    to    the 
reliability   shorthand    was  given    by   Rspicky[ Ref .2 ].      This 
paper    is   devoted    to    a    complementary    part    of   the    idea. 

Any  study   on    system  re-liability   always   requires   two 
steps;    one    is    the   description   of    the    system's   life   and   the 
other    is   the   derivation   of    its    survival    function. 

Under   the   assumption    of    constant    component    failure    rates 
this    paper    presents    a    way   of   obtaining    the    system's    reli- 
ability   which    requires   little   beyond    the    description    of   the 
system's   life. 

Section    II    contains   an   approach    to    the   convolution    of 
exponentially   distributed   random    variables.      Also    there    is   a 
presentation    of   a    computational    algorithm    for    the    convolu- 
tion  of   exponentially    distributed   random    variables. 


Appendix  A  gives  survival  functions  corresponding  to  sev- 
eral reliability  shorthand  notations  and  a  program  writer,  in 
Fortran  for  computations  from  shorthand  notations. 

Section  III  deals  with  the  reliability  of  redundant  sys- 
tems under  the  assumption  of  constant  component  failure 
rates. 

Appendix  B  consists  of  a  program  written  in  Fortran.  The 
program  supports  the  approach  of  Section  III.  There  is  a 
crude  Monte  Carlo  simulation  program  in  Appendix  D  which  is 
a  simulation  program  parallel  to  the  program  in  Appendix  B. 

There  is  another  program  in  Appeniix  C  writtsn  in  For- 
tran. This  program  uses  the  network  approach  to  systems 
described  in  Section  III. 

Appendix  E  summarizes  the  definitions  used  in  this  paper. 


II.   AN  APPROACH  TO  COMPUTING  CONVOLUTIONS  OF  EXPONENTIAL 

EIHP3H "VISTA  BI Eg  


This  section  introduces  a  general  algorithm  for  com- 
puting the  survival  function  of  any  convolution  of  exponen- 
tial random  variables. 

In  reliability,  the  teem  convolution  refers  to  the  summa- 
tion of  independent  random  lives.   la  order  to  have  simplic- 
ity in  specifying  convolutions,  tha  reliability  shorthand 
introduces  a  special  notation.   In  the  following  sections  we 
will  use  this  notation. 


A.  THE  SURVIVAL  FUNCTION  FOR  A  CONVOLUTION  OF  RANDOM 
VARIABLES 


Let  F1(t)  and  F2(t)  be  the  survival  functions  for  the 
random  variables  T1  and  T2  respectively.  Let  f1  (t)  and  f 2  (t) 
be  the  corresponding  densities.  Let  F  (t)  be  the  survival 
function  for  the  random  variable  T,  where  T=T1+T2. 

Then  the  likelihood  expression  for  F (t)  is 

t 

F(t)=F1  (t)  +/F2(t-s)  f1(s)  is 

0 


In    the   right    hand    side    of   equation,    F1  (•■-)     is   the    prob- 
ability   that   component   one    completes    the    mission, 

t 
J*F2(t-s)  f  1  (s)  ds   is   the   probability   that   at   some    time 

s    (0<s<t)    component   one   fails,    component   two  takes   its   place 
and    carries   ~he   system  to   the   end    of    the    mission    duration    t. 

In   order  to   illustrate  consider   some   applications. 
1 .      Example 

Reliability   Shorthand  Notation    :    EXP  {Xx }  +  EXP  (X2 } 
SYSTEM:      One   component   having   one   spire    with   a   dissimilar 
failure   rate.    If   the  active   component   fails,    the 
spare   will   replace    it   immediately. 

Here   the  life    for   the   system    is    T   =    T1    +   T2.       The    reli- 
ability   shorthand   notation    indicates   that    this   system    has   an 
exponential   life   with    failure   rate    Xi    followed   by    an    expo- 
nential  life    with   failure   rate   X2  . 

The    survival   function    for   the  active   component   is 

F,   (t)=   e"Xlt  ,    t>0. 


The    survival   function    for  the   spare   is 
F2  (t)  =   e'A2t  ,    t>0. 


The    survival   function    for   the   system    is 
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F(t)      =        R     (t)         ♦/     F?(t«S)      ff.     (3)      Is 
I  0         c 

F(t)    =   e-Alt+fe  -**(*-*)   XliXlS  ds      ,    t>0. 

If    we  complete  the   integration,    tha    result   is 


F(t)  =  X2/(X2-  Xi)     eXlt  *Xi/  (Xi-X2)     e"A2t    ,    t>0. 


Which  is  the  well  known  result. 

Another  way  to  establish  this  foraula  is  the  use  of  the 
moment  generating  function.  (Freund  and  Walpole[  Ref  .3  ]) 
2 .   Example 

Reliability  Shorthand  Notation  :  EXP{ X  \  +EXP  {  a } 
SYSTEM:   One  component  having  one  identical  spare.   If  the 
active  component  fails,  the  spare  will  replace  it 
immediately. 

The  formula  that  we  derived  above  becomes  meaningless, 
because  the  denominators  become  zero.  If  we  proceed  as 

before 

t 
F(t)  =  F1(t)  +J*F2(t-s)  f1  (s)  ds 

0 

-Xt    rt.A(t-sX  -Xs 
F(t)  ■  e   +  J  e  v   Xe        ds 


F(t)  =  (  1  ♦  \t    )  eAt     ,  t>3. 
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The  result  comes  out  as  expected  to  be  the  Erlang  (2,X  } 
survival  function. 


B.  RELIABILITY  SHORTHAND  NOTATION 

EXP{A}  +  EXP  {A}  + +EXP{A}   (n  identical 

exponential  lives) 

This  leads  to  the  Erlang  {r,  A.  }  survival  function 

_      n  x+ 

P(t)  =  I    (At)  /(i-1)!   e~Al  ,  t>3. 

i=l 

The   use    of   moment   generating    function   gives   the  result 
immediately. 


C.     THE    RELIABILITY    SHORTHAND    NOTATION 

EXP{Ai}  +  EXP{A2}  + *EXP{X    } 

This    is    the    expression   for  the   convolution   of    n   random 
variables    where    each    has   a    distinct    failure   rate. 

By   the   approach    used   in    Section    2.1,    adding   one    exponen- 
tial   life   at    once,    one  can    obtain   the    formula    for    the    surv: 

val    function 
n 

F(t)=   I     n.  a.  /  n    (Xl-x2)     e"Ar    ,  r^o. 

i=l    jjtl     J       j.^1 
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This    is    also    a    well  known   formula    and   can    be    obtained 
from    by    a   moment   generating    function. 

D.    THE    RELIABILITY    SHOP.TH&ND    NOTATION 

EXP{Ai]  ♦EXP(A2}  + +  EXP{A2}  (n4    terms) 

♦    EX?{A2}+EXP{A2}  + +  EXP{\2}  (rz    terms) 

+    EXP{Ak}  +  EXP[Ak}  + +EXP{Ak)  (nk    terms) 

k 
This   is    the  convolution    of   £      n.      exponential    random 

1-1       " 
variables   where  there   are  n^     identical   exponential   random 

variables  having  the   faililre  rate   X . . 

The   moment   generating   function   technique   is   not    useful   in 
this    situation,    since    there    is  a    hugs    amount   of   complexity 
involved.    This   section   deals    with   this    notation    using   the 
convolution   formula. 

1.      Be liability   Shorthand  Notation 

EXP{Ai}  +  EXP{Ai}+EX?  [Xz] 

We   know    the   survival   function    for    EXP  [\l } +SXP {A2}    and 
also    for  EXP{A  }+2XP{A}.    The   next   two    subsections   present 
different    ways   to    reach    the    survival    function    for    the    short- 
hand   notation   above. 
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(a)  . 

Let   T1,T2,T3    be   random    variables    distributed    as 
EXP  {Xi},    EX?{Ai},    EXP{A2}    respectively. 

Let   T1,=T1+T2.    This  random  variable    has   the   Erlang    {2,  A  } 
distribution. 

The  convolution    formula    for   T=T1»+T3    is; 


FT(t)=FT/(t)     + 


t 
f    _ 


FT    (t-s)    f^   (s)     ds         ,    t>0. 


-Ait. 


FT(t)  =  (1  +  X!t)e  "1W+ 


1 


0 


-A,t 


F    (t)=e      [(X22-2\1A2)/(A2-Ai)2         ♦   AiA2/( A  z-\  i  )     t} 
T 

*e'X2t    {     A12/(Xl-A2>  2    }     ,t>0. 
(b). 

Let   T1,T2,T3    be   random  variables    distributed   as 
EXP{Ai}f    EXP{Ai}#    EXP  (A  2}    respectively. 

Let   T2*=T2+T3.    This  random   variable    has   the   survival 


function 


-A,t 


FT    (t)    =X2/(A2-Ai)e  +A1/(A1~A2)e 

'2 


A2t 


,t>0. 


The   convolution    formula    for   T=T1+r2*    is 


FT(t)     =   FT<    (t)     ♦    f    FT      (t-s)     fT   ,(s)     ds        ,    t>0. 

I  i2  j      i1  i2 

0 
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FT(t)    =  A2/(A2-AX)    eXlt*X1/(X1-X2|     eXat*f  5XliJ^>(x2-X1J 


-    -Al  S         -A2  S.         a 

(e        -e       )     ds 


-Ait 


Bj  (t)     =e         {(A2  2-2  AiA2)/(A2  -\x  )  2    ♦  Ai  /(A2-Ai)     t} 
+eA2tU12   /(Ai-A2)  2    },t>0. 


Subsections    (a)    and    (b)     illustrate   that   the  convolution 
formula   gives   a    unique  result,    regardless   of   the    way   of 
choosing  the   prior   random  life. 

2.      Reliability   Shorthand  Notation 
EXP  (Ai)  *EXP  £Xi  }+EXP  [A  2}+EXP  {A2  } 

Let    T1,    T2,    T3,    T4    be  random   variables  distributed 
as 
EXP{Ai},    EXPfAi},    EXP{A2},    EXP(x2}    r?spectively . 

Then   from   the   derivation    in  Section    2,    T1'=T1+T2+T3    is    a 

random   variable   having  the    survival   function 

Fr,   (t)    =    (a   +a   t)eAlt^a2ie"Ai^t>0. 
Ui  11     12 

where    an  =  U22-2AiA2)  /  (A2  "  Ai)  2    ,    al2=  AjAy  (\2  "Xl  )     , 
a2i=A12/(A1"A2  )2    • 


The   convolution    formula    for  T=T1«+T4    is 

-        f- 

FT  (t)     =    FT/(t)+  JF    (t-s)fT,(s)ds    ,    t>0. 

FT(t)     =(aii+ai2t)eAlt*a21eX2t+|iX2(t(Hi^ra12+A1a123)e"Xl 

^  .  0 

+  A2  a2^  2     )    is 
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The  result    from   the   above   is 
F  T(+.)  =  (a^af*)  eAlt+  (4*^t|  e"Alt ,  t  >0 . 
where      a£=  (X?-3  X^Aj  )  /(A2  -  Ax)  3      ,    V2=  A:A2  2/  (A2  -  h  )  z      , 
a2,i=(A13-3A12A2)/(A1  -A2)3       ,    a^-Aa  2  A2/  (  Ar  A2  )  2       - 

It    is   important   to   noxe    that    the   number   of   exponential 
terms    equals   the   number  of    dissimilar    failure  rates   and  each 
exponential   term   has   a   polynomial   coafficient   with   the 
degree  of   the   polynomial   equal   to   the    number   of   identical 
random   variables   having  the    corresponding   failure    rate. 

The  next   section   deals   with   the   convolution   of    exponen- 
tially distributed   random   variables    asing   the   fact   illus- 
trated above. 

3.      Introduction   of   an    Algorithm    for   the  Convolution   of 
Exponential   Lives" 

From   Subsection   2,     we   can   infer    the   fom    of   the   sur- 
vival   function    which    we   sought  at   the    beginning    of 
Section   C. 

SHORTHAND    NOTATION    : 

EXP  {A  J+EXPCAi  }♦ +EXPCAJ}  (r.L     terns) 

+    EXP{A2}  +  EXP{A2}  + +EXPCA2}  (nfc    terras) 


♦    EXP{Ak}+EX?{Ak}  + +EXP{Akf  (nk    t^ras) 
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SURVIVAL    FUNCTION    : 

F(t)    =    A!(t)eA2t+    A2(t)    iA2t+...-f    A     (t)     eV   ,    t>0 

n  -i 

where   Ai  (t )  =a  n  +a  12  t+  a     tz+ ..  .. +a2n    t         , 

13  ,l2 

n2-i 
A  2  (t)  =a2i +a  22t+a  23t2+ ..  ..+a2n   t  , 

2 


K. 


a.       Example: 

SHORTHAND   NOTATION    : 
EXP{Xi}  +  EXP{X2  } 

*  EXP{X2) 

♦  EXP{A3  }+EXP{X3}  +  EXP{A3  } 


SURVIVAL    FUNCTION    : 

F(t)  =  (a1i+a12t)eAlt  +a21  sXlt  +  (a31  +  a  32  -  +  a33 t2)  eXz\  t>0 


17 


b.   An  Algorithm  to  Compute  the  Coefficients 

The  algorithm  represented  below  develops  the 
survival  function  by  adding  one  random  variable  in  each  run. 
As  an  example  ,    in  order  to    compute  the  survival  function 
for  the  convolution  of  ten  exponentially  distributed  random 
variables,  the  algorithm  is  supposed  to  run  ten  times. 

The  notation  used  in  the  algorithm  is: 

K     number  of  dissimilar  failure  rates 

X-j     ith  type  failure  rate 

X.  failure  rate  for  the  currently  entering  life 

'• 

a.  kth   coefficient   on  the    ith    polynomial 

n.  current   number    of    identical    lives   having 

the   ith   failure    rate 

nn.         number   of   random    variables   having   ith 
i 

failure   rate. 
Initial:  aik=0    »  V  J#k   where    j  =  1,2,..K      ,    k  =  1 , 2  , . .  .  r.j 

n.j=0    ,   V-      where   i=1,2,..K 
Input    :        Xj     ,  v-   where   i=1,2,...K 


nn. , V  i   where   i=1,2,...K 


The    first   run   is    :         n-|=1,       an=1 
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Algorithm    : 


n.=n.--M,    until   n,-  =nn,-    . 
]e     \  \       \ 

1.    Update  the  coefficients   :    a.   kfor    k=2,3,..,n, 

e  e 

ai  k=^iai      /  <m  — 1)     ,    where    m=n.    -j      for   j  =  0,1,.,n.  -1 
e        e   em-1  e  ne 


2.    Update  the   coefficiant    :    a. 
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3i    =ai  }*J  (a^Xi/cAe^)*  I     (j-1)!^  a-1j/(A2-tyJ} 


e        e    i^i 


n.>0 


j=2 
nf>l 


3.    Update  the  other  coefficients:      a.,     Vi*k 

where   i#i    , n . *Q   for    i=1,2,.,K,    k=1,2,.,n. 
e       l  I 


a.     =a.    A.  /(A.  -A.)    ,  vi  where   i*i    . 
ini      'Ve       'a     * 

ain    =(X1ain-ma1mH)/(XrV    '    vi    where   i#ie     and   V1 
e  e 

for      m=n.-j,     j=1 ,  2r . . ,  n..  -  1 


c.       Example: 

Reliability    Shorthand   Notation    : 
EX?{AX  }  +  EXP{a1}*EXP£x2}^SXP{A2} 

Let   T1rT2,T3,T4    be    random   variabiss   distributed   as 
EXP  {Ai},    EXPCAi),    EXP{A2},     EXP{A2  }    respectively. 
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We    would   like   to   derivs    a   formula    for   the   survival    func- 
tion  of   the  random   variable    T=TH-T2  +  r3  +  TU. 

The   use    of  algorithm: 
1 st    RON; 

n2  =1,    an  =1 

At    the   end   of   this    run   we   have   only    one    random    variable, 
which    is   distributed    EXP[AA}. 

2nd  SON; 


i      =1,    n=2 

e  -  1 


ai2 =  X1a11  then     a12  =    \x 

an  =a  ii  +  0  then     a  1X=1 

At   the   end   of   the    2.      run   we   have   the   survival    function 
for    the   T»=T1+T2,    where   T1    and   T2    are    identically    distri- 
buted   as   EXPCXi).    The    survival   function    is 


F(t)  =  (a11+a12t)e-X^ 


F(t)=  (1+ ALt)  e      ,    t>0 


Which    is   EBLANG (2, Xi} 


3  ^   RON; 


i   =2,    n   =2,    n    =1 
e 
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a21  =0>a11A1/(X1-X2  )+a       1!    X2/(\l'\2)z 
=Ui/(Ai-  h)  +  XiX2/(Aa-A2)  2 

=  Ai2/(A1-X2)2 
ai2  =a12A2/(A2  -Ai)         then        a12  =AXA2/  (A,  "  Ax ) 

an  =(A2aira12)/(A2  -Ax)       then    aai  =  (A23-2A!A2)  /  (  A2-Ax  )2 
Note   that,    here   the   coefficient    a  12    is    the   updated  one. 

At   the    end   of   the    3.      run    we    have   the    survival    function 

for    the   random   variable   T"-T1+T2+T3,    where    T1,T2,T3   are   as 

defined    before. 

-Ait        -A2t 
^(tJMaifa^tJe     +a2Je      ,    t>0 

where   a11=  <\z*-2^\2)  /  {\2  -\x  )  2  ,    ^\1x/(\2  "\x  )  > 

a2r=Ai2/(Ai  ~A2 ) 2  • 

Note    that   the   coefficients   are   identical   to   the    result    of 
Subsection     (b)  . 


4th    RUN; 


i   =2,    n,  =2,    n2=2 

e  l  2 


a22=A2a21  then        a22  «  \xz  \2  /( A1*A2  ) 

a21  =a21  +a11x1/(Ai  ~A2)*ai^2/(A  i~  A2)  2 
then        a21  =(Al3    -3X22x2  )  /  (^  -  A2>  3 
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ai2=ai2X2/(X2"X1)  --he*         a12  =A1X22/(A2-Xi  ) 


1i  =  tA2a1rai2)/<A2-X1)       the*      a1:  =U2  3"3x1A22)/(A2-A1) 


a  ;  = 


At  the  end  of  the  4.   run,  we  have  the  survival  function 
for  the  variable  T=T1 +T2+T3 +TU. 

FT(t)  =  (aljfa12-)iAlt+(a21fa22t)IA2t  ,t>0 
where    a12  =(A2  3"3A1A22)/(A2-AI)3    ,      a  12=  Xl A22/  (a2"  Ax  )  2    ' 

a2i  =  (A13-3A12A2)/(A1"A2)3    »      a22  =  *i  2  A2  /  ( Xi  "A  2  )  2    . 

Note   that   the  coafficisnts   are   identical  to   the   result   of 
Subsection    2. 
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III.       RELIABILITY    SHORTHAND    APPROACH    TO    SYSTEM    RELIABILITY 


This    section    deals   with  a    system    whose   components   have 
constant    failure   rates. 

Having   the   reliability   network   for   a    system   in    which   each 
component   has   an   exponential   life   and    knowing   the    probabili- 
ties   for   failures    (discussed   in   Appendix   S2)    makes   it    easy 
to   describe  the   system's   life. 

In    order   to    make   the   idea   clear,    we    will   go  through   some 
examples. 

A.     EXAMPLE 


2   spares    (Ai) 


SYSTEM    :    2    components    in    series   with    spares. 

DESCRIPTION    :    System   has    2   components    in    series.    Compo- 
nent   1    has    a   life   distributed  as   EXP{Xi}    and  there   are    2 
identical   spares.    Component    2   has   a    Life   distributed    as 
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EXP{A2}    and  there   is    one   p. on-identical    spare,    whose   life    is 
distributed  as  EXP{A3}. 
LIFE   : 

EXP{Ai+A2} 


/ 

3    *i/(Ai+X2) 

EXP(  A!+A2}  + 

1  h/lX^Xz  )  »     A2/(Ai  +  A2) 

/  \ 

EXP{A1  +  A2}+  EXPCXi+Aa  )  + 

2  A2/(Ai+A2) 

,ZERO 


EXPfAi  +A3  } 
7    Xi  /(Ai  *A3) 


A/(Ai  +A3  ) 

\ 

ZERO 


6    A3/(Ai<\\3) 

/ 

EXP^+As}^  /EXP{Ai+A3} 

5    Ai/(A!+A3)  9  Ai/(Xi   +A3) 

NEXP{Xln3  }\ 
\ 

10  A3/(A!+A3) 
^ZERO 

Explanation    for   the  derivation   of    system^   life: 

At  the  beginning  the  system  has  an  EXP{Ai+A2}  life.  The 
failure  of  the  system  is  by  the  failure  of  component  1  with 
a  probability  of  Ai/(Ai+A2)  or  by  tha  failure  of  component  2 
with    probability   of  A2  /  (Xi  +A2). 


2'4 


In  the  life  figure,  number  1  denotes  the  event  "failure 
of  component  1"  and  number  2  denotes  the  even*  "failure  of 
component   2". 

If  event  1  occurs,  component  1  is  replaced  by  one  of  the 
spares  and  system  again  functions  with  a  life  distributed  as 
EXP{A1+A2},  since  the  exponential  distribution  has  the  memo- 
ryless  property  and  component  2  still  has  the  same  failure 
rate. 

If  the  event  2  occurs,  component  2  is  replaced  by  its 
spare   and  the  system   has   a    life   distributed  as   EXP{Xi  +A3  }  . 

The  numbers  on  the  life  figure  correspond  to  the  transi- 
tions that  can  occur.  The  probability  on  each  arc  shews  the 
conditional  probability  of  the  transition.  As  an  example 
event  3  can  occur  with  probability  of  Ax  /(\l*X2),  given  that 
event    1    has   occured   before. 

The   distribution   ZERO   defined   by    Esary  [  Eef .  1  ]    and 
Repicky  C  Ref  .2]    (also      defined   in     Appendix   S3)         enters   when 
life    is    exhausted. 

For  convenience  of  description,  it  is  helpful  to  define 
the    concept      of    path   used   in     this   paper.      Path      denotes   the 
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sequence   of      events   in  the      system's    life    from      the   starting 
point    to   the   point    where   the   system   is   not    functioning. 
Examples; 

Events  1  and  3  are  a  path,  which  denotes  a  sequence  of 
lives  for  the  system.  In  this  case,  the  system  has  3  expo- 
nentially distributed  lives  EXP  {\x  ♦  A2 }  +  EXP  C^i  *A2  *  +sxp  (\i  *X2 } 
with    the    probability   of         \±    /(\1*X2)     Ai/(A!+A2>- 

Events  1,  1  and  8  form  another  path,  which  describes  a 
sequence  of  lives  for  the  system.  In  this  path,  the  sys-em 
life  is  EXP{  A  x +A2  } +SXP(Ai  +A2  )  +  EXp(Ai*  A3  }  +ZERO.  The  prob- 
ability  of    this    path    is 

[Ai/(Ai*A2][a2/(A!  +A2J[a3  /(Xi   *A3  ij. 

The  ZERO  distribution  contributes  zero  additional  life  to 
the  system,  so  we  can  omit  it.  Nevertheless,  we  can  not  omit 
its  probability  in  the  calculation  of  path  probability, 
since  the  event  numbered  3  has  a  probability  of  occurring. 
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B.     THE    SURVIVAL    FUNCTION    FOR    THE    SYSTEM 

In   this    section    we    will    deal      with    an   approach   to    obtain- 
ing   the    system's   survival   function. 
1 .       EXAMPLES 

SYSTEM:    Two    components   in    parallel. 


LIFE 


,Ai   /(X1*X2) EXP(x2} 


2XP{Ai  +  \2}     * 


\ 


Xz    /(Ai+Xz) SXPCA^ 


DESCRIPTION:    Component   1    has   EXP^}     life    and   component    2 
has    EXP(A2}    life. 

There  are  two  failure  events  that  can  be  allowed.  The 
first  is  that  component  1  fails  at  some  time  t  and  component 
2   carries   the    system    for   the  rest    of    the   time.    The   ether   one 
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is  that  component  2  fails  at  some  time  during  the  mission 
and  component  1  carries  the  system  for  the  rest  of  the  time. 
There  are  more  events  such  as  no  failures  during  the  mission 
duration,    which   are   taken  care   of   by   the    representation. 

Now   we   have   two    paths 


I°-°f    Path  l§JL2ht 


Life 


1  Pi  =Ai/(Ai+A2>  EXPUi+Aal+EXPftj   | 

2  P2=A2/(x1+A2)  EXP{Al  +X2}  +  EXP{Xl} 

Let   T  be  the   system's    time   to    failure    and   let    T1,T2,T3    be 
random   variables   exponentially   distributed      with    the    failure 
rates    \x*  \z,\2,  X2    respectively. 
The  n 

f  T1    ♦   T2  with   probability    p4 


m     _ 


( 


!fcT1    +    T3  with    probability    p. 


The   survival   function   oan   ba    writtan   as 

?(t)   =  pj   r,  (t)    +  ?2    f2  (t)      r  t>o 

where   p1  =  \2  /  (  X1+A2  )  '    P2  =  A2  /  (Ax  +  A2  )     *n3    F,  (t)  ,    F2  (t)       denote 
the    survival   functions   for    the   shorthand   notations 
EXP  (X1  +A2}  +  EXP{A2}  ,    EXP(Ai  *  A2  }+EX?CAi  }     respectively. 
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The  survival  function  for  the  convolution  of  two  exponen- 
tially distributed  random  variables  with  dissimilar  failure 
rates  is 

F(t)  =  W(A2-Ai)  eXlt+X1/(Ai-^  )  3A2x>0. 

If   we   do   the    substitutions  for   ?i  (t)     and   F2  (t)    as 
X\   s\i  +\z  ,  X2   -  Xz  andAi=Ai+A2,A2  =  X%     respectively, 
?1  (t)    becomes 

P1  (t)=  A2/(-Ai)     e(Al+X2)t+(A1+A2)/A1  eA2t    ,    t>0, 

and   F2  (t)    becomes 

-(Ai+A2)t  -A2t 

\   (*)  =A1/("A2)     3  +  (A1+A2)/A2e  ,    t>9. 

Then 

F(t)  =  X1/(A1  +  A2)     Fl  (t)     +  A2/  (Ai  +A2)     ?    (t) 

F(t)  =  eAl     +eXl     -iXl+A2)t    ,    t>0. 

The  result  gives  the  survival  function  that  is  wall  known 
for    this  system. 

Another    example    is    the   2    out    of   3    system. 
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RELIABILITY    NETWORK 


1 

2 

The    known  survival    function    is 

-Ui+A2)t  -(A^A3)t  -GW-Aalt      -(A!+A2+A3)t 


F(t)=e  ♦  «= 


♦  a 


-2s 


,    t>0. 


LIFE; 


,EXP£x2+A3  * 


A2  /(A!  +  A2  +  A3) 


EXP  {A2  +A2+Aa  }  ♦  —  \2/<Ai  ^A^Aa  )  — EX?{xx  +A3) 


\/  (Xi+Xz  *A3) 


•EXP{Ai  +  A2} 


'EX?  (  Ai  ♦  A2+A2  }  +  EXP  {A2  +  A3  }    with    probability  \x/  (A^A^) 


T=  i  EXP{  A!  +  A2+A3  } +EXP  {Ai+A3}     with    probability  x2/  (A  i%%) 


EX?  fXi  +  XZ  +  X3  }  +EXP  CA1  +  A2  }    with    probability  \/  (xf-fc+x} 
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The    survival   function    is 


F(t)=p3  ?!    (t)+p2F2    ft)*P3   ?3  (t)     ,    t>0. 


where 

-(X!+X2+X3)t  -(X2+X3)t 

Fi  (t)  =  (\2+X^  /  (-XL)    e  ♦  (XX<X2  fX3)  /h  « 

-(X!+X2+X3)t  -CXi+X3)t 

F2  (♦-)  =  (X^  /  (-X2)     e  ♦  (X/X,+X^  A2  •      ■ 

-(Xi+X2+X3)t  -(Xi+X2)t 

F3  (+-)  =  (Xi+Xs)  /  ("X3)     e  +  (X/Xz+X^  A 3  e 

If  we  do  the  necessary  cancellations,  we  can  get, 

p(t)se-  (A2+x3)t+e-(x1+x3)t  +e-(Xx+x2)t  _2e-(x1+x2+x3)tj  t>0 


as  desired. 

2-   General  Procedure 

Having  the  algorithm  presented  in  Section  II,   we  can 
treat  more  complicated  systems  under  similar  assumptions. 

The  procedure  is 

i.   Set  up  the  reliability  network  for  the  system. 

ii.  According  to  this  network,  set  up  the  system  life. 

iii. Using  the  proper  reliability  shorthand  formulas  for 
the  related  convolutions  of  exponential  lives,  set  up  the 
survival  function. 
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IV.        SUMMARY 


The  reliability    shorthand   is    an   easy    way   to    describe   a 
system's   life,    but    it    is    difficult    to    implement    ccmputatio- 
naly    since    there   is   considerable   complexity   in   handling   con- 
volut  ions. 

The   algorithm   presented    in   this   paper   gives   5oa>e   relief 
from    this   difficulty.    However,    the   accuracy    in   obtained    from 
this    algorithm   is   very   much    related    to    the    differences    in 
the    failure   rates. 

Another  aspect  in  the  algorithm  i  =  that  distributions  are 
convolved  one  at  a  time  and  this  requires  very  accurate  run- 
ning   conditions    in    the  case    of   a      complicated    system. 

It    is    believed    that   it    is   possible    to    derive   another 
algorithm   which    is    more    powerful  than    the   one    introduced 
here.       Instead   of   adding    one   distribution    at    a    *-.ime,    one   can 
try    to   add    several   at    a    tine. 
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APPENDIX       A 


A-1 


This  section  contains  the  survival  functions  for  several 
shorthand  notations  which  were  derived  by  the  use  of  the 
approach   described    in    Section  II. 


A. 1.1    Shorthand    Notation    :    5XP{A}. 


Survival  Function:    F(t)=   e  ,    t>0. 


A.  1.2.1    Shorthand    Notation:    EX?(Ai}  +    EXP{A2} 


Survival  Function; 


F(t)=  A2/(A2-Ai)     eXl    +Ai/(A1-A2)     eA2t  ,t>3 


A.  1.2.2    Shorthand   Notation    :    EXP { X } +EXP  {  \  } . 


Survival   Function: 


F(t)=  (1+  A  t)  e"At       ,    t>0 


A.  1.3.1    Shorthand    Notation    :    EXP  {\x }  +EXP  {\z  }  +EX?(x3  } 


Survival   Function: 


F<t)=a21eAlt*a21iA2t*a3i    **A3t       ,    t>0. 
where   axl  =a2X3/(a2  -\i>   h3'X^     '    *21  =  AiA3/ ( Ai  -A2  )  (A3  -  A2  )     , 
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a3  1   =  X1X2/  (  Xl  ~  A3  )     (X2~X3)» 


A.  1.3.2    EXP(X   }+EXP{  X }  +  EXP(  X  } 


Survival  Function: 


F(t)  =  (UX  t+1/2X2t2)  9"At       ,    t>0 


A.  1.3.3    EXP{Xi  }+EXP(X2} +EXP{X2  } 


Survival  Function: 


F(t)=ane"Xlt+(a      *a      t)e"A2t    ,    t>0 

21        22 


where   au=  X22/  (  \2'\l  )  2  ,    %=  (Xx  2-2X2X2)  /  <Xi  -X2)  2 

2  2~  Xi  Xy    V  X  i    a  2  1       * 

A.  1.4.1    EXP{X!}>EXP{X2  }  *  EXP {  A3  }+EXP{X-  } 
Survival   Function: 

F(t)=axleXl  ♦a2ieX2t+a31"^3t*a»reX,»t    ,    t>0. 

where      a.,  ■     H     Xj/  n     (X.-X.)        v  i=1,2,3,4 
11       J7i         j/M       J 

A.  1.4.2    EXP(  X  }+EXP[  X}*-  EX?{X   }  +  EXP{  X  } 
Survival   Function: 


F(t)  =  (1+\t*1/2  X  2t2  +  1/6   A3t3)iAt     ,    t>0 


A.  1.4.3    EXP{Xi}  +  SXP{X2  }  ♦  EXP(X2}+E!CP{X2  } 
Survival  Function: 
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-Ait  -\2t 

F(t)=a11«      +(a21+a22t+^3t2)s      ,    ^>0 


where   an=  \2  3/(x2-A!  )  3      ,    \{*^-\x 

a22=A1A2(Af2x2)/(X1-A2)  2      *    \f  Xx  A2  2/2  (Xj  -Xj  )  . 
A. 1.4.4    EXPCXi  }+EXP{A1}*EXP(A2}+EXP[A2  } 

Survival  Function: 

-A2t  -A2t 

F(t)  =  (a1*a12t)e      *(a21+^2t)s      ,    t>0. 

where      au=  (  A23-3\  2\1 )  /(x2  -  a  l )  3,    a12=  A^2/  (A2  -  Ax)  2  r 

a2i=(Ai3-3A2A2)  /(Ai  -A2)  3r    a22=A2Ai2/(A1-A2  )2. 

A.  1.4.5    EXP  (A!  }+EXP{A2}  ♦  EXP{  A  3} +EX?  {A3  } 

Survival  Function: 

-Ait         -A2t  -A3t 

P(t)=a   e      +a   e      +  (a   +  a    t) s        ,    t>0. 

11  21  31       32 

where   d^  =  ^A3  /(  \2-  Ai )  (A3-Ai  )  2,    ay  =  A^/fAi"^  )  (  A3  "A2  ) 

a3=  AxA2/(  A1-A3)   <A2  -A3)     +A1A2A3[  1  /  (  Ax  -,\2  )  (Ai  -\3>  2 

-1/(Aj  -A2)  (A2  -A3  )  2     ] 


a    «  Ai^As/  (  Xl  "  A3)   (  A  2"  A  3 )  • 
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A- 2 

This  section   introduces  a  Fortran   prograa   using  the 
algorithm  described  in  Section  II. 

A  PROGRAM  FOR  THE  ALGORITHM  ONE  AT  A  TIME: 


c 
C 
C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 


THIS    IS    A    PROGRAM    TO    COM 

WHICH    HAS    THE    RELIAB 

EXP  (L1) *. 

WHERE  THERE  IS  NO  RESTR 


VARIABLES: 

I.TH  TYPE  FAI 

ON  TH 

AMOUNT  OF  LI 7 

FAI 

AUXILAURY  ARR 

FOR 
THE  ARRAY  FOR 


A(I,J) 

NI(I) 

NINIT(I) 

MI) 

REAL  L(20),A(20.10) 
INTEGER  NINIT  (20)  ,  NI  ( 


PUTE  THE  RELIABILITY  OF  A  SYSTEM 

ILITY  SHORTHAND  NOTATION 

.♦EXP  (LN) 

ICTION  FOR  THE  FAILURE  RATES. 


LORE  RATE,  J.TH  COEFFICIENT 
E  POLIN0M. 

ES  HAVING  THE  I.TH  TYPE 
LURE  RATE. 

AY  FOR  THE  NI  (I) .  INPUT 
THE  PROGRAM. 
THE  FAILURE  RATES 
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20) 


GET    INPUT 


READJ5.499)     T 
WRITE(6,498)     T     . 
READ(5,500)     K 
WRITE  (6.1500)     K 
READ  (5,501)      (NINIT  (I) 
WRITE(6,1501)      (NINIT( 

CD 

KK  =  1 

A  (KR,  1)  =1  . 

NI  (KK)=1 

JJ=1 

IF(JJ.EQ.KK)     30 
A  (KKf  1)  =A  (KK,1) 
IF  (JJ.EQ.K) 


GO    TO 


TO    8 


JJ=JJ+1 

GO    TO    4 
IF(KK.EQ.K)     GO 
KK=KK+1 
GO    TO    2 
CONTINUE 


IE=1 

IF  (NINIT(IE)  .  EQ.NI  (I 

NI  (IE)=NI  (IE)  +1 


,L(I)  ,1=1,  K) 
I)  ,L(I)  ,1=1,  K) 

MPUTE    COEFFICIENTS    ONE    ALL 
DISSIMILAR 


TO    6 
*L(JJ)  /(L(JJ)  -L  (KK)  ) 

9 


BEGIN    TO    ADD    ONE    AT    A    TIME 


E))     GO    TO    99 


UPDATE       IS 
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c 
c 

c 


c 
c 
c 


J  =  0 
12  NNNN=NI(IE)-J 

A  (IE,NNNN)  =L  (IE)  *A 
IF  (NNNN.EQ.2)     GO    T 
J  =  j+1 
GO    TO    12 


(IE.NNNN-1)  /FLOAT  (NNNN-1) 
0    20 


20    SUM=0. 


1/K 


E)     GO 


>  * 


DO    101    1= 

IF(I.EQ.I 

SUM=SUM+A 

IP(SIfl)  .  IT'.  2)     GO*    T 

FACT=1. 

NKK=NI  (I) 

DO 

FAC 


TO    101 

LjI!/(L(I)  -L  (IE)) 


102 
101 
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sum=sum+a  (i 
:on 


102    II 

T=FACT 
,11) *F 
TINUE 


=2,NKK 

<II-1) 

ACT*L  (IE)  /  (L  (I)  -L  (IS)  )  **II 


CONTINUE 

A  (IE,1)  =SUM+A(IE,  1) 


UPDATE      I.NE.IE 


1=1 
IF  (I 
A  (I, 


EQ.IS)     GO 


NI(I))=A(I, 

F  (Nlfl)  .EQ.  1)     GO 


TO    25 
NI(I) 


\j-f 


*L(IE)/(L(IE)-L(I) 
TO    26 


2a 

NNI  =  NI  (I)  -J 
A(I,NNI)=(L 

*                                         / 

IF(J.EQ.  (NI 
J=J  +  1 

GO    TO    2U 

26 

IF(I.EQ.K)     GO    TO 
1  =  1+1 
GO    TO    21 

99 

CONTINUE 

IF(IE.EQ.K)     GO    TO 
IE=IE+1 
GO    TO    32 

c 

c 

c 

IE) 


*A  (I,NNI)  -FLOAT  (NNI)  *A  (I,NNI+1)  ) 

E)-L(IT) 

1)  )     GO    TO    26 


32 


35 


41 

40 


35 


CONTINUE 

PRO=0. 

DO    40    1=1, K 

SUM=0. 

NKK=NI(I) 

DO    4  1    J=1,NKK 
SUM=SUM+A  (I,  J)  *T*< 
CONTINUE 

PRO=PRO+SUM*EXP  (-L  (I) 

CONTINUE 


CALCULATE  THE  PROBABILITY 


J-1) 
*T) 


PRINT    OUT    THE    RESULTS 


17 


498 
499 
500 
501 


WRIT 
DO  1 
NNK= 
WRIT 
WRIT 
WRIT 
STOP 
FORM 
FORM 
FORM 
FORM 


E  (6,600) 
7    1=1, K 

NI(I) 
E  (6, 1600) 
E  (6,603) 
E(6,601 


L  (I) 

T,  pr5 


,NI 
J)  , 


AT 
AT 
AT 
AT 


I1) 
J=1,NNK) 

MSIOXj'THE    INPUT    IS:  '  ,  //1  5X  ,  «TIMS=  '  ,?  10.  4) 

F10.4) 

3X,I3 

3X,I3,F10.5) 
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600  FORMAT(»1  ■  ,  10X,  'FAILURE    RATS  AMOUNT      COEFFICIENTS 
*■  ,/10X,40  ('-')) 

601  POBHAT("0' ,10X, 'PROBA BILITY    (    TIME    >    ',F10.5,')  =    'r 
*F10.  7) 

6  03    FORMAT  ('0  •  ,  35X,  3(5  (E1  2.  5,2X)  ,/>  ) 

1500  FORMATJ'O' , 14X,  'K         = »  .  3X, 12 ,7X , •  (THE    NUMBER    OF 
♦DISSIMILAR    FAILURE    RATES)  ») 

1501  FORMAT('0' .15X, '#    OF     R.V.  FAILURE    RATS', /15X, 

^10." 

13) 


125  ('-')  ,//20M 8 X.I 3,5X,F10.5,//)  ) 
1600    FOEMAT(»0',ldx, F10.5,  5X,I~ 


END 

EXAMPLE  : 

INPUT  FOR  THE  PROGRAM 


K=    6 

5  6.5 

5  1.5 

5  2.5 

5  3.5 

5  5.0 

5  9.0 


OUTPUT    FROM    THE    PROGRAM 


THE    INPUT    IS: 
TIME=  8.0000 


K         =  6  (THE    NUMBER    OF    DISSIMILAR    LAMDA"    S) 

NO.    OF    R.V.  FAILURE    RATE 


5 

6.5 

5 

1.  5 

5 

2.5 

5 

3.5 

5 

5.0 

5 

9.0 

FAILURE    RATS 

AMOUNT 

COEFFICIENTS 

6.5  5 

0.15713E+G8  0.84332S+07         0.1353QE+07 

0.20038E+06  0.92701E+0'4 
1.5                             5 

0.27364E+09  -0  .  89949E+03         0.11817E  +  08 

-0.73813E+06  0 . 1 8586E+05 
2.5                             5 

-0.18295E+11  3.82312E+10       -0.17601E+10 

0.17341E+09  -0.11979E+08 
3.5                            5 

0.18841E+11  0.89164E+10         0.25648E+10 

0.25077E+09  0.334  06E+03 
5.0                            5 

-0.83480E+09  -0.52324E+09       -0.16061E+09 

-0.21713E+08  -0.22945E+07 
9.0                             5 

-0.39700E+05  -0.15168E+06       -0.32323E+06 

-0.33727E+06  -0.27860E-01 

PROBABILITY     (    TIME    >  9.0    )=       0.6935042 
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APPENDIX  B 

This  section  contains  a.  program   described  in  Section  III 
to  compute  the  reliability  of  a  system. 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS  PROGRAM  CALCULATES   THE  RELIABILITY  FUNCTION 


?(T)  = 
P1       *       F1  (T 

P2      *      F2  (T 
P3      *      F3  (T 


EXP  (L1r1)  +EXP 
'EXP  (L2,1)  +  EXP 
•EXP  (L3,1)  +SXP 


L1,2)  + +EXP  (L1  ,MV 

L2,2)  + +  EXP  (L2,H2! 

L3,2)  + +  SX?  (L3VM3 


PI       *      FI  (T)  (  =  EXP  (LI,1)  +  SXP  (LI,  2)  + +EXP  (LI ,  MI)  ) 


PN      *      FN  (T)  (=EXP  (LN,  1)  +EXP  (LNf  2)  + +  EXP(LN,MN)) 


PROBA 
TIME 
A  (I,K) 


RATE 


ARRAY 
ARRAY 
NUMBE 
TOTAL 
PROBA 
TIME 
THE  C 
I.T 
AND    K 


1 


1 

6 

10 
1 


c 
c 

c 


REAL    P( 

CALL    RE 

SUM=0. 
DO    1 
SUM= 

IF  (ABS( 

PROBA=0 
DO  2 
WRIT 
CALL 
PROB 

WRITE  (6 

STOP 

WRITS  (6 

STOP 

FORMAT ( 
*F10.7) 
09    FORMAT( 
98    FORMAT( 
*    EQUAL 

END 


50) 

AD1  (IK,P 

1=1, IK 

SUM  +  P(I) 

SUM-1.)  . 
-0 

1=1, IK 
E(6,1009 

ONECON( 
A=PROBA+ 
,601)     T, 


FOR    ALL   LAMDAS. 

FOR    ALL    PROBABILITIES. 
R    OF    EX?    IN    EACH    ROW. 

NUMBER    OF    ROWS. 
BILIT?    OF    SYSTEM    SURVIVAL    AT    TIME    T. 

©EFFICIENT    FOR    THE    CURRENT    PATH, 

H    TYPE     FAILURE 

.TH    COEFFICIENT    ON    THE    THIS    POLINOMIAL 


,T) 


GT.  1.0E-5)     GO    TO    199 
)    1,1, P(D 

pro, if 

PRO*P  (  I) 

PRO 


99 
01 


,198)     SUM 

•  1  •  ,10X,  'PROBABILITY    (    TIME    >    ',F10.5,')  =    ', 

12,'. •  ,3X.'P(',I2  •) 
•THE  SUM  OF  THE  PRO! 
/10X, ' SUM    =    ',F10.7) 


f0'  ,10X,I2,«.  •  ,3X,  •  P('  ,12,  •)  ='  ,F10.  7) 
M«,10X,  'THE    SUM    OF    THE    PROBABILITIES    IS    NOT 
TO    1.0«, 


GET    INP0T 


SUBROUTINE 
REAL  P(50) 
READ  (5,499) 


READ1  (IK,P  ,T) 
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3 

499 

500 

1499 

1500 


C 
C 
C 


WRIT 
READ 
WRIT 
READ 
RETU 
FORM 
FORM 
FORM 
FORM 
FORM 
END 


S(6,1499)     T 
50  0)     IK 


ft 


,1500)     IK 
(5,3)      (P(I)  ,1=1, IK) 
RN 

AT(7F10.7) 
AT  '5X,F10.4) 
AT(5X,I5) 

AT('1  ■  ,  10X,  'TIME=  •    F10.4} 
aT(v0v,10Zv «IK      =  «3X,I3) 


SUBROUTINE 
REAL    A(20,10), 
INTEGER    NI  (20) 


C 
C 

c 


1 

52 
1501 


ONECON(PR0,T) 

INTEGER'    NI  (20f  ,  NINIT  (20) 

CALL  READ (K,L,NINIT) 

CALL  ONEDIS  (Kf ArNIrL) 

CALL  ONEATA  (K,A,NI,L,  NINIT) 

CALL  CALPRO  (K,A,NI,L,  PRO, T) 

RETURN 

END 


SUBROUTINE  READ  (K,  Lr  N  INIT) 
REAL  L(20) 
INTEGER  NINIT  (20) 
READ  (5,1)  K 
WRITE  (6,1501)  K 

READ(5,52)  (L  (I)  ,  NINI  T  (I)  ,1  =  1,  S) 
RETURN 

FORMAT(5X,I5) 
FORMAT(5X,F10.5,I5) 
FOPMAT(10X, »K    =',I2f5X, •(*  OF 
*  RATES.) •) 
END 


DISSIMILAR  FAILURE 


C 
C 

c 


SUBR 
REAL 
INTE 
KK=1 
A  (KK 
NI  (K 
JJ  =  1 
IF  (J 
A(KK 
IF  (J 
JJ=J 
GO  T 
IF  (K 
KK=K 
GO  T 
RETU 
END 


OUTINE    ONEDIS(K,A,NI,  L) 

A(20  ,10)  -L  (20) 
GER    NI(20J 

,1)  =1. 
K)  =  1 

J.EQ.KK)     GO    TO    6 

,1)  =A  (KK,  1)  *L(JJ)  /(L(JJ)  -L(KK)  ) 

J.2Q.  K)     GO    TO    9 

J  +  1 

o  4 

K.EQ.  K    )     GO    TO    3 
K  +  1 
O    2 

RN 


BEGIN    TO    ADD    ONE    AT    A    TIME 


SUBROUTINE    ONEATA(K,  A,NI,  L, NINIT) 
REAL    L(20),A(20  ,10! 

er  nin: 


(IT  (2 


i% 


1(20) 


32 


C 
C 
C 


INTEG 

IE=1 

IF  (NINIT  (IE)  .  EO.  NI  (IE))     GO 

NI(IE)=NI(IE)  >1 


TO    9  9 


UPDATE      IE 


40 


J  =  0 
12  NN  =  NI{IS)-J 

A  (IE.NN)  =L  (IE)  *A(IE.NN-1)  /FLOAT  (NN-1) 
IFfNN.  EQ.2)     GO   TO    20 
J=J+1 
GO   TO    12 
20    CO NT I  NO E 

DO    101    1=1, K 

IF(I.EQ.IE)     GO    TO    101 

AUE,lf=A  (IE.1)  ♦  A(I,1)*L(I)/(L(I)-L(IE)) 

IF(NIfI).LT.  2)     GO    TO    101 

FACT=1. 

NKK  =  NI(I) 

DO    102    II=2,NKK 
?ACT=FACI(II-1) 
A{IE,1)=A  (IE,  1)+A(I,II)  *FACT*L  (IE)  /(L(I)  -L  (IS)  )  **II 
102  CONTINUE 

101  CONTINUE 

C 
C  UPDATE      I. HE. IS 


C 


1=1 
21 


IF(I.SQ.IE)     GO    TO    26 
A(I,NI(I)  )=AJI,NI(I)  )  *L(IE)/(L(IE)  -L  (I)  ) 
IP(NIfl).EQ.I)     GO    TO    26 


\j=t 

24  NNI  =  NI(I)-J 


ii 


A(I,NNI)  =  (L(IE)  *A(I  .NNIj  -FLOAT  (NNI)  *A  (I,NNI+1)  ) 
/  I  Li  I X  Si)  •"  L  I XI  ) 
IF  (J.SQ.  (NI  (I)  -  1)  )     GO    TO    26 


J=J+1 
GO   TO    2U 

26    IF(I.EQ.K)     GO    TO    32 

1  =  1+1 

GO  TO  21 
99    CONTINUE 

IF(IE.SQ.K)     GO    TO    35 

IS=IE+1 

GO  TO  32 
35    FETURN 

END 


SUBROUTINE    CAL?RO(K,A  ,NI,  L,PRO,T) 

REAL    L(20)  , A  (20  ,10) 

INTEGER    NI(20) 

PRO=0. 

DO    40    1=1, K 

SUM=0. 

NKK=NI(I) 

DO    41    J=1,  NKK 
SUM=SUM+A  (I,  J)  *T*<  J-1) 
41  CONTINUE 

TTT=-L(I)  *T 
~RO=?RO+SUM*SXP  (TTT) 
40    CONTINUE 

WRITE  (6,1010) 

DO    100  0    1=1,  K 
NNN=NI  (I) 

WRITE  (6, 1001)     L(I)  ,NI(I)  ,  (A  (I,  J)  ,J=1,NNN) 
1000  CONTINUE 

RETURN 
10  01    FOPMAT(10X,F10.5,5X,  I3,5X,  .2(5  (E1  2.  5  ,5X)  ,/30X)  ) 
1010    FORMATJ15X, 'LAMDA    ', 5 X, ■ NI ■ ,3 X, • COEFFICIENTS ', /1 UX , 

135C-')) 

END 
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INPUT    POI     THE    PROGRAM     : 


TIMS 

15. 

#ROWS 

5 

.2 

.2 

1.ST 

3 

.  1 

3 

.8 

2 

.4 

1 

2.  ND 

3 

.2 

2 

.5 

4 

.6 

1 

3.  RD 

3 

.  1 

3 

.5 

4 

.3 

2 

4.  TH 

2 

.  1 

3 

.  2 

3 

5.TH 

3 

.  1 

2 

.2 

4 

.4 

3 

.2  .2 


OUTPUT  FROM  THE  PROGRAM 


TIME  =         15.0000 
IK       =  5 

1.         P  (1)  =   0.2 

K         =3  (#    OF    DISSIMILAR    FAILURE    RATES.) 

LAMDA       HI       COEFFICIENTS 


0.10         3         0.11294E+01    0.6S342S-01    0.87075E-02 
0.80         2         0.  18742E-G1    0.23324E-02 
0.40         1       -0.14815E+00 

2.  P  ( 2)  =   0  . 2 

K         =  3  (#    OF    DISSIMILAR    FAILURE    RATES.) 

LAMDA      NI      COEFFICIENTS 

0.20         2      -0.  25077E  +  02   0.23  148E  +  01 

0.50         4      -0.  13017E+03    0.21333E+02    -0 . 22222E-01    0.55556E-01 

0.60         1         0.15625E  +  03 

3.  ?  (3)  =    0.2 

K  =3  (#    OF    DISSIMILAR    FAILURE    RATES.) 

LAMDA       NI    COEFFICIENTS 


0.10         3         0.75531E+01    -0.54932E+00       0.27466E-01 
0.50         4      -0.89945E  +  01    -0.99536E+00    -0.42847E-01    -0.73242S-03 
0.30         2         0.24414E  +  01    -0.14648E+01 
4.  P  (    4)  =    0.  2 

K         =2  ^#    OF    DISSIMILAR    FAILURE    RATES.) 

LAMDA       NI    COEFFICIENTS 


0.10         3       0.32000E+02    -0.16000E  +  01       0.40000E-01 
0.20         3    -0.31000E+02    -0.14000E+01    -0.20000E-01 
5.         P  (    5)  =    0. 2 
K         =3  i#    OF    DISSIMILAR    FAILURE    RATES.) 

LAMDA       NI    COEFFICIENTS 


0.10         2   -0.  15170E+03  0.37926E  +  01 

0.20         4      0.14400E+03  0.12800E+02    0.32000E+00    0.10667S-01 

0.40         3      0.87037E+01  0.51851E+00    0.38389E-02 

PROBABILITY     (    TIME  >    15.0    )=       0.9111188 
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APPENDIX  C 
INTRODUCTION 

This  section   consists  of  a   computer  program   to  compute 
system  reliabilities  as  described  in  Section  III. 

Again  for  simplicity,  we  will  go  thrucgh  an  example. 

The  structure  of  the  system  is 

spare  ( Xx ) 


Components  1,  2  and  3  have  lives  exponentially  distributed 
with  failure  rates  \i  ,\2  ,  A3  respectively.  Also  we  have  a 
spare  for  component  1  which  is  identical  to  component  1. 

Using  the  shorthand  approach,   the  system  life   would  b<= 
determined  as  follows 
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Xi   /(Xi+A2+A3) EXP{A2+A3} 

.EXP  {\+^+A3}+-A2  /^i+A2+A^ SXPCAi+As) 

A^CV^+A^  A3/(V1+A2+A3) ZERO 

EXP  {A1+A2+A3}  +-A3/  (A^+y  -ZERO 

A2/C^*A2+A3)  Ai/C^+^J-SXP^+Ag} 

^EXP{A  +  ^       *f 

\ 

As/fA^J-ZSRO 

Seme   definitions   are   necessary   before   describing  the   pro- 
gram 

NODE:  Each  node  represents  an  exponentially  distrubuted 
random  variable  with  a  certain  failure  rate.  The  number  for 
any  node  can  be  chosen  arbitrarily  but  can  not  be  used  more 
than   once. 

ARC  :  Each  arc  originates  at  a  node  and  leads  to  another 
node.  Only  one  arc  can  terminate  at  a  given  node.  Arc  num- 
bers can  be  chosen  arbitrarily.  If  n  is  the  total  number  of 
nodes  in  a  system  life  then  there  will  be  n- 1  arcs  in  this 
system. 

With   these   definitions,    we  can   represent   a   sys-em    life    in 
the    following   way 
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3ACK  POINTER  LIST  (IP B) :  Each  node  has  a  back  pointer.  A 
back  pointer  is  an  arc  number  which  shows  which  arc  connects 
the  node  to  the  tree.  If  the  pointer  is  zero,  then  the 
related  node  is  the  root  of  the  tree. 


NODE  NO 
1 
2 
3 
a 
5 
6 
7 
8 
9 


BACK    POINTER    IPB(I) 
0 
1 
2 
3 

5 
6 

7 


(Here  node  1  is  the  root  of  the  tree.) 


NODE  CODE  LIST:  As  we  mentioned,  each  node  represents  an 
exponential  lifetime,  for  simplicity  we  can  use  some  integer 
code  numbers  for  each  failure  rate. 
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Code   No.  Relatsd  Failure    Rate-L(U 

0  (distribution    ZERO) 

1  Xi  +X2  +  X3 

Xi        x3 


3  Xx    +     X 


2 


Here  only  the  code  number  0  is  not  arbitrary  and  0  can  be 
used  for  the  ZERO  distribution.  The  others  can  be  picked  ou* 
arbitrarily. 

ARC    CODE  LIST:         This    list    is    similar    to    the   node    code    list. 
Arc   code   numbers   represent    probabilities. 
Code    No.  Prob  abilz^ty^PA  (I) 

1  \l/CX1*X2*X3) 

2  X2  /{Xi+X2+X3  ) 

3  X3/ai+X2+X3) 

*  Xi/(Xi+X3) 

5  X3/(Xi  +  X3) 


ARC  ORIGIN  LIST  :  Each  arc  has  an  origin  node  and  a  termi- 
nal node.  In  the  program  we  need  to  use  only  the  origin 
list. 


U6 


Arc   No. 
1 


Q^isio.  il2l§  IQ  (I) 


4 
5 
6 
7 
8 


LAST  POINT  NODE  LIST  (LP):  This  list  indicates  the  nodes  on 
the  end  of  each  path.  There  is  no  necessary  order  in  the 
list. 

LAST  POINT  DEPTH  LIST  (IPD)  :  This  list:  indicates  the  number 
of   arcs    from   last    node  to   the   root    of    the   tree. 


I    L.  ?.  Node  List  LP  (I) 
1  5 


2 
3 
4 
5 
6 


6 
7 
8 
9 
4 


k-  P.  DeRlk  List  IPD (I) 
2 
2 
2 
2 
2 
1 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS    PROGRAM    COMPUTES    THE    RELIABILITY    OF    A    SYSTEM 
WHICH    HAS    A    SURVIVAL    FUNCTION    LIKE    BELOW: 

P1       *      F1  (T)  (  =  EXP  (L1  ,1)  +  EXP  (L1.2)  + ♦  EXP(L1.M1)) 

♦  P2       *       F2(T)   (=EXP(L2r  1)  «-EX?(L2,2)  + ♦EX?(L2,M2) 

+  P3       *       F1(T      =EXP(L3,  U  *EXP(L3,2)  *■ +SXP(L3,K3 


♦    PI       *      FI  (T)  (  =  EXP  (LI,1)  +EXP  (LI,  2)  + +2XP  (LI  , MI) ) 

+ 

4-  .       .  .       . 


+    PN      *      FN  (T)  (  =  EXP  (LN,1)  +EXP  (LNf  2)  + +£XP(LN,MN)) 


IPB  (I) 

IPA(K) 

IO(K) 

ICL(I) 

EMI 

I    ° 

N 

M 

NLAST 

NDIF 

MDIF 

NPATH() 

HPATH() 

IM 

IP 

PROBA 

PRO 

P 

LCODE() 

LL() 

IK 

NIO 
NINIT() 
A  (,) 
T 


VARIABLES 
BACK 
PROBA 
ORIGI 
FAILU 
DIFFE 
LAST 
LAST 
FAILU 
NUMBE 
NUMBE 
NUMBE 
NUMBE 
NUMBE 
NODE 
ARC  L 
THE  N 
THE  D 
THE  C 
PROE. 
THE  C 
FAILU 
FAILU 
#  OF 
NUMBE 
AUXIL 
COEFF 
TIME 


IN  PRO 
POINT E 

BiLirr 

N  MODE 
RE  RAT 
RENT  P 
POINT 
POINT 
RE  RAT 
R  OF  N 

OF 

OF 

OF 

OF 


A 
L 
D 
D 

LIST  I 
1ST  IN 
UMBER 
EPTH  0 
URRENT 

of  ra 

URRENT 
RE  RAT 
RE  RAT 
DIS.  ? 
R  OF  L 
AURY  A 
ICIEST 


GRAM: 
R    OF 

CODE 

OF    K 
E    COD 
ROBAB 
NODE 
NODE 
E    LIS 
ODES 
PCS    I 
AST    P 
IFFER 
IFFER 
NAP 

A  PA 
OF  TH 
F    THE 

VALJ 
E    COR 

PATH 
E  COD 
E  LIS 
AILUR 
IFES 
RRAY 
S    FOR 


THE  I 
OF  r 
.TH  A 
E  OF 

iLiry 

LIST 

DEPTH 

T 

IN  IS 

N  TRE 

OINT 

ENT  F 

ENT  P 

ARTIE 

RTica 

E  CUR 

co  a  9 

E  OF 
RENT 

PR  OB 
ES  IN 
T  IN 
E  RAT 
RELAT 
FOR  N 

CORK 


. TH  NODE 
HE  K.TH  ARE 
RE 

THE  I.TH  NUDE 
LIST 


LIST 


OD 

IL 

OB 

LA 

AR 

EN 

NT 

HE 

AT 

31 

CO 

OR 

S 

D 

NT  PATH 


ES 
OR 

AB 
R 

p 

T 

L 
S 
H 

LI 
RR 
RE 
IN 
TO 


S  RATES 

ILITIES 

PATH 

ATH 

LAST  POINT  NODE 

AST  POINT  NODE 

URVIVAL  FUNCTION 

SURVIVAL  AT  T 

TV 

ENT  PATH 
NT  PATH 

CURRENT  PATH 
THE  EAILORE  RATE 


REAL 

INTE 

♦NPAT 

PROB 

CALL 

WRIT 

D 

I 

I 

NPA 

DO  9 

NPAT 

MPAT 

IM=N 

i  CONT 


MAIN  PROGRAM 

,  (20)  ,?A  (2  0)  ,  LLM1)  ,  PRO,  P  .  PROB  A,  T 
:a  IPB  (100)  ,ICL  100)  ,10  (13  0)  ,LPJ50)  ,IPD(! 
"11)  ,  MP  ATH  (10)  r  SI  (11)  ,1,  I*,  I  PA  (100) 
=0.0 

IE  AD (LP ,L,  PA, IPB  ,ICL , IO# IPD, NLAST ,T, I ? A) 
[6,101  3) 
1  1=1, NLAST 
=LP(I) 


:H  (1  1)  ,MPATH  (10)  , 

SA=0.0 

,    READ(LP,L, PA, IPB ,1 

:E  (6,1013) 
JO    1    1=1, NLAST 
IM=LP(I) 
IP=IPD  (I) 
TH  (1)=IM 
9    J=1,I? 

TH  (J  +  l)=IO  (IPB(IM)  ) 
TH  (J)  =I?B(IM) 
"PATH  (J+1) 
INUE 


U3 


CALL    PASS /PA, L, ICL, NPATH, M PATH, IP, LL,NI, IK,  ?,IPA) 
KRITS(6,1990      1,1, P 
CALL    ONECON(PRO,T,  LL,NI,IK) 
PROBA=PROSA+P*PR0 
1  CONTINUE 

WRITE(6,601)     T,PROBA 
STOP 
501    FORMAT (»01  ,10X,  'PROBABILITY    (    TIME    >    '  ,F10.5  ,')  =    ', 
*F10.7)  l  '    '  ' 

1313    FOPMATM0X,  'BEGIN    TO    CALCULATION'     /1  1X,  20  ( '  -  • )  ) 
1990    FORHAT(«0*  ,10X,I2,'.'  ,5X,  «P(»,I2,M  =  '  ,F10. 7) 
END 
C 

c  get  iNpar 

c 

SUBROUTINE  READ  (LP,L,  PA,IPB, ICL , 10, IPD, NLAST, T , IPA) 

REAL  L(20),PA(20) 

INTEGER  IPB(100)  ,ICL(100)  ,10  (1)0)  ,  LP  (50)  ,  IPD  (5  0)  , 

C  READ   TIME  AND  #  OF  NODE 

C 

READ  (5,  1)  T,N 
C 

C  HEAD  BACK  POINTSRLISI, FAILURE 

C  PATE  CODE  LIST 

C 

READ  (5,3)   (IPB(I)  ,  ICL(I)  ,I  =  1,N) 
M  =  N-  1 
C 

C  • READ  ARC  ORIGIN  LIST  AND  ARC 

C  PROBABILITY  CODE  LIST 

C 

READ  (5,4)   (IO(K)  ,1  PA  (K)  ,K=1,  M) 
C 

C  READ  #  OF  PATH  AND  #  OF  DIFFERENT 

C  FAILURE  RATE  AND  #  OF  DIFFERENT 

C  ARC  PROBABILITIES 

READ  (5,5)  NLAST, ttOIF,MDIF 
C 

C  READ  LAST  POINT  LI3T  (LP)  AND 

C  LAST  POINT  DEPTH  LIST  (IPD) 

C 

READ  (5,6)   (L P  (I)  ,  IPD  (I)  ,1=1,  NLAST) 
C 

C  READ  FAILURE  RATES 

C 

READ  (5,7)   (L  (I)  ,1=  1,NDIF) 

READ  (5,7)   (PA(I)  ,I=1,MDIF) 
RETURN 
1  FORMAT(5X,F10.3,I5) 

3  FOPMAT(5X,10I5] 

4  FORMAT(5X,10I5J 

5  F0RMAT(5X,3I5) 

6  FOPMAT(5X,10I5) 

7  FORMAT(5X,5F10.  5) 
END 

C 

C         THIS  SUBROUTINE  COMPUTES  NECESSARY  LISTS  TO 

C   CALCULATE  THE  SURVIVAL  FUNCTION  OF  A  CURRENT  PATH 

r 

C 

SUBROUTINE  PASS  (PA , L, ICL, NPATH, MP ATH , IP , LL , NI , IK, F, 
*IPA) 
REAL  PA  (20)  ,LLM1)  ,L(20) 

INTEGER  NPATH  (1  1)  ,  «PA  TH  (1  0)  ,  LCD  DS  ( 10)  ,  NI  ( 1  0)  ,  ICL  ( 1  00)  , 
*I?A  (100) 
P=1. 
DO  1  1=1,1? 
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1  P=P*PA(IPAr«PATH(I)  )  ) 

IPP=IP-M 

11=1 

I?  (ICL  (NPATH  (1)  )  .SQ.O)     11  =  2 

DO    2    I=II,IPP 

2  LCODE(I)  =ICL  (NPATH  (I)) 

IK=0 

DO    3    1  =  11  ,I?P 

IP(LCODEfl)  .EQ.  0)     GO    TO    3 

IK=IK+1 

LL  (IK>  =L(LCDDE(I)  ) 

LLL=1 
JJ=I+1 
IF(JJ.GT.IPP)     10    TO    13 

DO    4    J=JJ,IPP 

IF(LCODE(I)  .  NE.LCODE  (J)  )     GO    TO    4 

LLL=LLL+1 

LCODE(J)  =0 
4  CONTINUE 


13  NI(IK)=LLL 

1(1)  =0 

CONTINUE 


LCODE 


RETURN 

END    C 
C  THIS    SUBROUTINE    CONTROLS    THE 

C      COMPUTING    OF    THE    SURVIVAL    FUNCTION    OF    THE    CURRENT    PATH 
C 

SUBROUTINE    ONECONfPRO  ,T,L, NINIT , K) 

REAL-  A  (11  ,10)  ,L  (1T) 

INTEGER  NI(11)  ,  NINIT(  11) 

CALL  ONEDIS  K, A,NI,L) 

CALL  ONEATA(K,A,NI,L, NINIT) 

CALL  CALPRO (K,A,NI,L, PROf I) 

RETURN 

END 
C 

C  THIS    ROUTINE    COMPUTES    THE    COEFFICIENTS 

C  AS    EACH    IS    DISSIMILAR 

C 

SUBROUTINE    ONEDIS  (  K,  A  ,NI,  L) 

REAL    A(11  ,10)  ,L  (11) 

INTEGER    NI(20) 

KK  =  1 
2    A  (KK,  1)  =1. 

HI  (KK)  =  1 

JJ  =  1 
4    IF  (JJ.EQ.KK)     GO    TO    6 

A(KK,  1)  =A  (KK,  1)  *L(JJ)  /  (L(  J  J)  -L  (KK)  ) 
6    IF  (JJ.EQ.  K)     GO    TO    9 

JJ=JJ  +  1 

GO    TO    4 
9    IF  (KK.EQ.K    )     GO    TO    3 

KK=KK+1 

GO    TO    2 
8    RETURN 

END 
C 

C  BEGIN    TO    ADD    ONE    AT    A    TIME 

C 

SUBROUTINE    ONEATA(K,  A,NI,  L, NINIT) 

REAL  L  ( 11)  ,A (11  ,10 

INTEGER  NINIT  (1  1)  ,NI  (  11) 

IE  =  1 
32  IF  (NINIT(IE)  .EQ.NI  (IE))  GO  TO  99 

NI  (IE)  =  NI  (IE)  +1 

C  UPDATE   IE 

C 
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c 


J=0 
12  NN=NI(IE)-J 

AjfIE,NN)  =L  (IE)  *A(IE.NN-1)/FL0AT  (HN-1) 
IF(NN.EQ.2)     GO    TO    20 
J=J+1 
GO    TO    12 
20    CONTINUE 

DO    101    1=1.  K 

IFJI.EQ.IEf    GO    TO    101 

MIE,lf=A  (IE.1)  +A(I,1)*L(I)/(L(I)-L(IE)) 


IP(NI(I)  .  LT.  2)     GO    TO    101 
FACT=1. 


NKK  =  NI  (I) 

DO    102    11=2, NKK 

*  FACT=FACT<II-1) 

A  (IE,  1)  =  A  (IE,  1)  +  A(I,II)  *FACT*L(IE)/(L(I|  -L  (IE))  **II 
102  CONTINUE 

101  CONTINUE 

C 
C  UPDATE      I.NE.IE 


1  =  1 
21    IF(I.EQ.IE)     GO    TO    25 

A(I,NI(I)  )=A(I,NI(I)  )  *L(IE)/(L(IE)-L(I)) 
IF  (Nim  .EQ.  1)     GO    TO    26 

24  NNI  =  NI  (I)  -J 

A  (I, NNI)  =  (L(IS)  *A  (I,  NNI)  -FLOAT  (NNI)  *A  (I,  NNI+1)  ) 


*  /(L(IE)-L(I)) 

IFfJ.SQ.  (NI(I)  -1)  )     GO    TO    26 
J=J+1 

GO    TO    24 
26    IF(I.EQ.K)     GO    TO   32 
1  =  1+1 
GO    TO    21 
99    CONTINUE 

IF(IE.EQ.K)     GO    TO    35 
IE=IE+1 
GO    TO    3  2 
35    RETURN 
END 
C 

C         CALCULATES    THE    PROBABILITY    FOR    T3E    CURRENT    PATH 
C 

SUBROUTINE    CALPRO (K, A , NI, L,PRO,T) 

REAL    L(11)  ,  AM1  ,10) 

INTEGER    NI(1T) 

PRO=0. 

DO    40    1=1, K 

SUM=0. 

NKK=NI(I) 

DO    41    J=1, NKK 
SUM=SUM+A  (I,  J)  *T*<  J-1) 
41  CONTINUE 

TTT=-L(I)  *T 
PRO=PRO+S0M*SXP  (TTT) 
40    CONTINUE 

WRITE  (6,1011) 
DO  1000  1=1, K 
NNN=NI  {I) 

WRITE  (6,1010)     L  (I)  ,NI   (I)  ,  (A  (I,  J)  ,J=1,NNN) 
1000    CONTINUE 

RETURN 
10  10    FOF.MAT(,0,,10X,F10.7/  5X,I  2  ,  5X,  2  (5  (E1  2.  5  ,2X)  ./2  5X)  ) 
10  11    FORMAT  (15X,«  LAM  DA'  ,7X  ,'NI'  ,5X,  •  COEFFICIENT  *  ,/10X, 

135('-M) 

END 
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INPUT 

DATA 

FOB   TH 

E    PROGRAM: 

T.N 
N(5.  1 

2. 

13 

0 

1 

1           1 

2 

1 

3 

3 

4 

1 

NO.  6 

5 

1 

6           2 

7 

3 

8 

3 

9 

2 

NO.  11 

10 

3 

11           2 

12 

2 

AS.  1 

1 

1 

1           2 

2 

1 

2 

2 

3 

1 

AR.  6 

3 

2 

4           3 

5 

1 

5 

2 

6 

1 

AH.  11 

6 

2 

7           3 

NL,  UN 
LPIPD 

6 

3 

3 

3 

3 

9           3 

10 

3 

11 

3 

12 

3 

LP6PD 

13 

3 

LAMDA 

2. 

0 

.5 

1.5 

PA'  S 

.25 

.75 

1.0 

OUTPUT    FROM    THE    PROGRAM    : 


BEGIN    TO    CALCULATION 

1  .  P  (    1)=    0.0625000 

LAMDA  NI    COEFFICIENTS 

1.5  2      -0.80000E+02        0. 24000E+02 

2.0  2         0.81000E+02         0.  18000E  +  02 

2.  P(    2)=    0.0  46  875  0 

LAMDA  NI    COEFFICIENTS 


1.5 
2.0 
3. 

LAMDA 


1.5 
2.0 

5. 
LAMDA 


-0. 51200E+03 
0. 51300E+03 
?  (   3)=    0. 1406250 


0. 96000E+02 
0.  16200E+03 


!Fi 


SI   COEFFICIENTS 


0.5  2      -0.63578E-06        0.  11852E+01 

2.0  3         0.10000E+01         0.81481E+00 

4.  P  (    4)  =    0.0468750 

LAMDA  NI    COEFFICIENTS 


2  -0.51200E+03 

3  0. 51300E+03 


NI    C 


P  (   5)=    0.  1406250 

FICIENTS 


0. 96000E+02 
0.  16200E  +  03 


OEF] 


0.5  2      -0.63578S-06        0.  11852E+01 

2.0  3         0.  10000E  +  01         0.81481E+00 
6.  P  [   6)=    0.5625000 

LAMDA  NT    COEFFICIENTS 


0.5  2      0.59259S+00 

2.0  2      0.40741E+00 

PROBABILITY     (    TIME 


0.  18000E  +  02 


0.22222E+00 


0.  18000E+02 


0.22222E+00 


0.38889E  +  00 
0.22222E  +  00 
2.0    )=      0.0173515 


52 


APPENDIX    D 

This  section  gives  tha  program  for  the  simulation  of  a 
system's  life  mentioned  in  Section  III.  The  program  uses  a 
crude   Monte   Carlo   simulation  procedure. 


C 
C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 

c 


c 
c 

c 


+ 
+ 
+ 

+ 
+ 


THIS    PROGRAM    SIMULATES    A    SYSTEM    HAVING 
THE    RELIABILITY    FUNCTION: 


?   T)  = 
P1       *      F1  (T 
P2      *      F2(T 
P3      *      F3(T 


=  EXP 

[LI,  1  )+EXP 

(L1,2 

=  EXP 

L2, 1J+EXP 

(L2,2 

=  EXP 

L3,  1     +  EXP 

(L3,2 

♦ +SXP  (L1  ,M1 

+ +  EXP    L2,M2 

♦ +  EXP    L3,  M3 


PI       *       FI(T)   (  =  EXP(LI,1)+EXP  (LI, 2)  + +  EXP(LI,MI)) 


?N      *       FN(T)   (=EXP(LN,  1  )  +  EXP(LN,  2)  ♦ +  EX?  (LN,  MN)  ) 


|| ! 

MT(    { 
N 

PROB 
TIME 


ARRAY    FOR    ALL   LAMDAS. 

ARRAY    FOR    ALL    PROBABILITIES. 

NUMBER    OF    EX?    IN    EACH    ROW. 

FIST    NUMBER    OF    LAMDAS     IN    THIS    ROW. 

TOTAL    NUMBER    OF    ROWS. 

PROBABILITY    OF    SYSTEM    SURVIVAL    AT    TIME 

TIME 


REAL    L(500)  ,P(50)  r?R0B 

INTEGER    MI(50)  ,MT(50)   ,  MIDO,  MIUP  ,1  ,  J,  N 

IX=U56378 

CALL    READ  (L,N, TIME, MI ,MT, P) 

CHECK    FOR    SUM    OF    P'S    EQUAL    TO    1.0 


150 


TOT=0. 

DO    150    1=1, N 
TOT=TOT+P  (I) 
IF  (A3S(T0T-1. 


STOP 


)  .GT.  1.E-5) 
CALL  CdNTRO (L,N,TIME, MI ,MT ,P, IX. PROB) 
CALL  OUTPUT (L#N, TIME, MI, MT, P, PROB) 
STOP 
END 


SUBROUTINE  READ 


SUBROUTINE  RE  AD jL, N,  P IME, MI, MT,  P) 
REAL  L(500)  ,P  (50) 
INTEGER  MI  (50)  , MT(50) 
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C      GET    THE    TIME 
C 

READ  (5,505)     TIME 
C 

C         GET    THE    NUMBER    OF    ROWS 
C 

READ (5,501)     N 
C 

C      GET    ALL    P«S 
C 

READ  (5,504)      (P(I),I=1,N) 
C 

C      GET    ALL    LAMDAS    IN    THE    ORDER    OF    RD»     BY    ROW 
C 

MIDO=1 

MIUP=0 

DO    100    J=1,N 

IF(J.NE.1)     MIDO=MIUP+1 

MT  (J)  =  MIDO 

READ  (5,  520)     MI  (J) 

MIUP=MIUP+MI  (J) 

READ(5,503)      (L  (I)  ,  I  =  MIDO,  MIUP) 
100    CONTINUE 

RETURN 
501    FORMAT (5X, 15) 

503  FORMAT(5F10.5) 

504  FORMAT(5F10.7 

505  FORMAT(5X,F10.3) 
520    FORMAT(5X,I5) 

END 
C 

C  SUBROUTINE    CONTRO 

C 

SUBROUTINE    CONTRO (L,N ,TIME, MI , MT, P , IX, PRO) 

REAL    L(500)  ,P(50) 

INTEGER    MI  (50) , MT(50) 

SUM=0. 

DO    1    1=1. N 
NN  =  1000000*P  (I) 

CALL    SIMULA(NN,L,MI(I)  ,MT(I)  ,X, IX, TIME) 
1  SUM=SUM+X 

PRO=SUM/1 000000 . 

RETURN 

END 
C 

C         SUBROUTINE    FOR    SIMULATION 
C 

SUBROUTINE    SIMULA(NN,  L,M,MTT,Xf  IX, TIME) 

REAL    L(500)  ,RN(50) 

INTEGER    M,I,J,IX 

MMT=MTT  +  M-1 

DO    1  1    1=1, NN 

TEST=0. 

CALL    LEXPN  (IX, RN,M ,16307,0) 

JJ=0 

DO    111    J=MTT,MMT 
JJ  =  JJ+1 
111  TEST=TEST  +  RN  (JJ)/L(J) 

11  IF  (TEST. GE. TIME)     X  =  1.0+X 

RETURN 
END 
C 

C  SUBROUTINE    OUTPUT 

C 

SUBROUTINE    OUTPUT (L,  N ,TIM E, MI, MT , ?, PRO) 
REAL    L(500)  ,P  (50) 
INTEGER    MI  (50)  , MT(50) 
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c 

C      PRINT    OUT    ALL    THE    SYSTEM 
C 

WRITE  (6,600) 

MIUP  =  MI  (1) 


WRITE  (6.601)     P(1)  ,  (L  (I)  ,I=1,MI0P) 
DO    102    1=2, N 


+  MI(I)-1 

2)     P(I)  ,  (L(J)  ,J=MIDO,MIU?) 
102    CONTINUE 

WRITE  (6  ,603)     TIME,  PRO 
RETURN 
C 
C      FOPMAT    STATEMENTS 


C 


6C0    FORMAT(» 1" ,5X,' THE    SYSTEM    IS    :•) 

6C1    FORMAT^O' ,5X,' F(T)  =•     F10.7    •     *     '  ,  ■  EXP         L    =    ',5 

1  (F10.4,3X)  ,5  (//2(k,5m0.  4,3X)  )  ) 
6C2    PORHAT(»(P,9X    «+    »,Ft0.7    •    *    *,4EXP         L   =    ',5(F1 

10.4,3X)  .5  (//2OX,5(F10.4,3X)  )  ) 
603    FORMATMO*  ,5X' PROBABILITY    OF    SYSTEM    SURVIVAL    AT 


10.4,3X)  .5  (//26x,5(F10.4,3x)  ) 
03    FORMATMO'  ,5X'  PROBABILITY    0 
*T=*,F10.5,f    IS     ',F10. 7) 


END 
INPUT    FOR    THE    PROGRAM 


TIME 

15. 

#ROW 

5 

.2 

.2 

.2 

.2 

.2 

1.  ST 

6 

.  1 

.8 

.4 

.1 

.  1 

.8 

2.  ND 

7 

.2 

.2 

.5 

.5 

.5 

.5 

.5 

3.  FD 

9 

.  1 

.5 

.3 

.1 

.  1 

.5 

.5 

.5 

.3 

4.  TH 

6 

.  1 

.2 

.1 

.1 

.2 

.2 

5.  TH 

9 

.  1 

.2 

.4 

.1 

.2 

.2 

.2 

.4 

.4 

OUTPUT    FROM    THE    PROGRAM    : 

THE    SYSTEM    IS    : 

F(T)=    0.2    *    EXP         L    =    0.1  0.8         0.  "4  0.1         0.1 

0.8 

+       0.2    *    EX?         L    =    0.2  0.2         0.5         0.6         0.5 

0.5  0.5 

♦  0.2    *    EXP         L    =    0.1  0.5         0.3  0.1         0.1 

0.5  0.5         0.5  0.3 

♦  0.2    *    EX?         L    =    0.1  0.2         0.1  0.1         0.2 

0    2 

+       0.2    *    2X?         L    =    0*.  1  0.2         0.4  0.1         0.2 

0.2  0.2         0.4         0.4 
PROBABILITY    OF    SYSTEM    SURVIVAL    AT    T=       15.00000    IS       0.911241 
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APPENDIX  E 

This  section  reviews  some  notions  which  are  found  in  the 
references  for  this  paper. 

E. 1  Redundant  systems  with  exponentially  lived  compo- 
nents: 

In  reliability  analysis,  the  term  system  is  used  to 
describe  a  set  of  components  organized  to  perform  some  mis- 
sion. A  system  is  redundant  if,  in  some  fashion,  some  of  the 
components  involved  act  as  back  up  for  other  components. 

A  rcugh  definition  might  be  that  a  system  is  not  redun- 
tann  if  the  failure  of  any  one  of  its  components  causes  the 
failure  of  the  system,  and  is  redundant  if  on=  or  more  of 
its  components  can  fail  without  causing  the  system  to  fail. 
Thus  redundant  systems  have  the  property  that  they  can  suf- 
fer damage  through  the  failure  of  some  of  their  components 
and  still  survive.  (ESARY'  B  ef .  1  ]) 

E. 2  First  failure  in  a  set  of  exponentially  lived  compo- 
nents : 
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We  have  n  components,  9 ach  independent  from  the  others, 
we  want  to  compute  the  probability  that  the  j.th  component 
fails  first. 


P(j.th  component  fails  first)  =P  (T  <T.  ,\/i  *i*j) 

J 


-A.S 


oJ 


P(T.j  >r.,vi,H*jir  =s)  ,\eJ  d< 


=  f  Cn     p,  (s)    )  ]f.  (s>    is 

\°°  -A-jS       -A.s 

[   n    e      ]x.e J     ds 

n  r°°       n  jfciV 

■    V(    I    i.)         [    J  X.  ]  e  ds 


0 


J        1;1       1       o; 

n 


i=l 


i 


E. 3  Degeneracy  at  zero  (Zero  Distribution): 

Let  ZERO   be  the   name  for  the   distribution  of   a  random 
variable  that  is  degenerate  at  zero. 

If  P(To=0)=1,    then  we  say  that  ro   has  the  distribution 
ZERO,  or 

T  ^  ZERO  . 

o 


57 


LIST    DF    REFERENCES 


1.  Esary,    J.    D.     ,    Unpublished   Course    notes   and    Problem    sets 
for   OA    4302,    Naval    Postgraduate   School, 

Monterey,    CA.939U0,    1982 

2.  Repicky,    John    J.,    Jr.  <    An   In  troduct  ion    to    a   Rili^bili^v 
Shorthand.    H-    S.    Thesis,    Naval   Postgraduate   ScEool, 
HonYefeyT   CA. 93940,    1981. 

3.  Freund,    J.    E.    and    Waipoie,    R.    E.    , 

Maremat ical    Statistics,     Third    edition,    1980. 


58 


INITIAL    DISTRIBUTION    LIST 


No. Copies 

1.  Defense    Technical    Information  Center  2 
Cameron    Station 

Alexandria,    Virginia   223  14 

2.  Library,   Code   0142  2 
Naval   Postgraduate    School 

Monterey  California   9394  0. 

3.  Departaant   Chairman,    Code   55  1 
Naval  Postgraduate   School 

Monterey  California   9394  0. 

4.  Professor   J.    D.    Esary,    55Ey  1 
Naval  Postgraduate   School 

Monterey  California  9394  0. 

5.  Professor      T.    Jayachandr an,    53Jy  1 
Naval   Postgraduate   School 

Monterey   California   93940. 

6.  Lt.    J.    G.    Sadan    Gursel  1 
Inonu  Cad.    No.    10 

Kesar.   -    Edirne   /   TURKEY 

7.  ortaiogu  Teknik    Universitesi  1 
Yen   Eylem   Arastirmasi   Bolumu 

Ankara/TURKEY 

8.  Dz.    K.    K.    iigi  2 
Bakanliklar,    Ankara  /   IURKEY 

9.  Lt.    Deniz   CORA  1 
Ovak    Sitesi   27.    Blok    Da.     1 

Yeni   Levent,    Istanbul   /    TURKEY 

10. Deniz   Harp   Okulu   K.    Iigi  1 

Heybeliada,    Istanbul  /TURKEY 


59 


Thesis 

G868 

c.l 


1983H* 

Gursel 

Some  computer 
algorithms  to  imple- 
ment a  reliability 
shorthand. 


25  SEP  84 


2  79  U8 


Thesis  I983HI 

G868         Gursel 

c.l         Some  computer 

algorithms  to  imple- 
ment a  reliability 
shorthand. 


thesG868 

Some  computer  algorithms  to  implement  a 


3  2768  001  03695  7 

DUDLEY  KNOX  LIBRARY 


Brvawi 

■' ^'"!---.'. '  i  ■''■"'■  ■■' 

OuKSm 

-'■',''''■:'*!■''.''?*', 

!i '    ' ' 

''■l-vkii^  9 

Ml 

'  .,/>i)v  QOOJjK  ..'»*''' 

a»sli 

mm  wtoeS ' R ; 

M 

yV,".v.''V'-".'.u*v,V 

HUffiMM 


' '  '  ■  ' "  » '' ' •■ '  v-' ' :  ' SHS S&S 1  BUS  3 


WtW 


^™R8&^-(!^v":> 


JflsSflfl 


vvBv 


'•■'  Ivfi 


:flfflfflOa»M»ag«^ 


-^:;;;ilillt'' 


nm 


